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A recent experiment [Nadj-Perge et al, Science 346, 602 (2014)] provides evidence for Majorana 
zero modes in iron (Fe) chains on the superconducting Pb(llO) surface. Here, we study this sys¬ 
tem by scanning tunneling microscopy using superconducting tips. This high-resolution technique 
resolves a rich subgap structure, including zero-energy excitations in some chains. We compare the 
symmetry properties of the data under voltage reversal against theoretical expectations and provide 
evidence that the putative Majorana signature overlaps with a previously unresolved low-energy 
resonance. Interpreting the data within a Majorana framework suggests that the topological gap is 
significantly smaller than previously believed. Aided by model calculations, we also analyze higher- 
energy features of the subgap spectrum and their relation to high-bias peaks which we associate 
with the Fe d-bands. 


Building on advances in nanofabrication [1], engineer¬ 
ing topological phases by proximity in superconducting 
hybrid structures has come within reach of current ex¬ 
periments. A major motivation for realizing such phases 
are their non-abelian Majorana quasiparticles [2-4], and 
their subsequent applications. The underlying topo¬ 
logical superconducting phases can be realized in one¬ 
dimensional (Id) helical liquids contacted by conven¬ 
tional 5-wave superconductors [5-9]. Among the most 
promising platforms studied in experiment are semi¬ 
conductor nanowires [10-15], edges of two-dimensional 
topological insulators [16, 17], and chains of magnetic 
adatoms [18, 19]. While the proximity coupling to a su¬ 
perconductor is needed to induce a gap protecting the 
topological phase, it also has more subtle consequences. 
Magnetic interactions mediated by the superconductor 
can stabilize magnetic order in the Id system [20-23]. 
Conversely, the spin structure may affect the supercon¬ 
ductor. This is particularly apparent for adatom chains, 
where a band of subgap Shiba states [24-27] may strongly 
modify the low-energy properties of the system [8, 28-32] 
and possibly induce trivial zero-energy features at the 
chain end [33]. At strong coupling, the Id states bleed 
substantially into the superconductor, reducing the effec¬ 
tive coherence length at low energies [34]. 

Nadj-Perge et al [18] recently provided intriguing evi¬ 
dence for Majorana states in Fe chains on Pb(llO). Here, 
we present data on the same system employing scanning 
tunneling microscopy/spectroscopy (STM/STS) with su¬ 
perconducting tips (see also [18]). We show that the use 
of superconducting tips not only provides enhanced res¬ 
olution of the subgap structure, but also allows for addi¬ 
tional consistency checks on the interpretation of the data 
in terms of Majorana quasiparticles. Our observations in¬ 
dicate a subgap spectrum comprising a flat Shiba band 
and strongly dispersing Fe states. An interpretation in 
terms of Majorana states suggests that the induced gap 
is considerably smaller than previously believed. 

We carried out the experiments in a Specs JT-STM 
at a temperature of 1.1 K. Cycles of sputtering and an¬ 


nealing of a Pb(llO) single crystal (Tc = 7.2 K) resulted 
in an atomically flat and clean surface. We employed 
Pb-covered superconducting tips (see Ref. [35] for the 
preparation procedure), which provide a resolution be¬ 
yond the Fermi-Dirac limit [36-38] (in our measurements: 

70 /ieV). Fe chains were prepared by e-beam evapo¬ 
ration from an iron rod (99.99% purity) onto the clean 
surface at room temperature, similar to Ref. [18]. With¬ 
out further annealing, we obtained chain lengths of up 
to 10 nm (measured between chain end and interven¬ 
ing Fe cluster). Single adatoms and dimers were pre¬ 
pared by e-beam evaporation onto the cold sample in 
the STM (T < 10 K) with a density of 350 adatoms 
per 100 X 100 nm^. The differential conductance dl/dV 
as a function of sample bias was recorded using stan¬ 
dard lock-in technique at 912 Hz (subgap spectra: bias 
modulation Rmod = 15 setpoint V = b mV, I = 

250 pA; large-bias spectra: Vnod = 2 mVrms, R = 2 V, 
7 = 850 pA). 

The d//dV-signal measures a convolution of the den¬ 
sities of states (DOS) of tip and sample when the tun¬ 
neling rate is slower than the quasiparticle relaxation of 
the subgap states [39]. This is the case for all measure¬ 
ments presented in this paper, as the differential conduc¬ 
tance varies linearly with junction transparency. STS 
with a superconducting tip shifts all sample peaks by 
TAtip. The coherence peaks then appear as resonances 
at eV = ±(Atip + Asampie)- Subgap states of the sub¬ 
strate of energy 5 appear as d//dV resonances at volt¬ 
ages ± (Atip + e) due to tunneling between the occupied 
(unoccupied) DOS singularity of the tip and the unoc¬ 
cupied (occupied) subgap state. At biases below Atip, 
resonances occur at ± (Atip — e) which result from ther¬ 
mally excited quasiparticles in substrate or tip. At our 
experimental temperature, this effect is limited to small 
energies 5. Accordingly, zero-energy Majorana states are 
signaled by resonances at eV = TAtip [18], and we thus 
pay particular attention to biases near this threshold. 

Figure 1 shows subgap d//dV-spectra, recorded at the 
termination of six independent Fe chains. The chains 
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Figure 1. dJ/dV-spectra recorded at the end of six different 
chains (small arrows mark position on the chain). Chains in 
(a-d) are terminated by a small cluster, in (e,f) they have a 
sharp cut-off. The energy resolution in (e,f) is reduced (width 
of BCS resonance: ~ 330 /iV), due to non bulk-like super¬ 
conductivity of the tip. Atip in mV: (a) 1.36, (b) 1.42, (c,d) 
1.38, (e,f) 1.24. Chain lengths measured between chain end 
and cluster onset in nm: (a) 13.9, (b) 9.5, (c) 6.2, (d) 6.0, (e) 
7.7, (f) 4.0. Scale bars correspond to 4 nm. For spectra of 
(b-f) at negative bias see Supplementary Material [40]. 


in Fig. l(a,b,d) exhibit a clear zero-energy signature at a 
bias of +Atip, which has been interpreted as a fingerprint 
of a Majorana bound state [18] (for the determination of 
the value of Atip, see Supplementary Material [40]). The 
zero-energy feature is accompanied by two resonances at 
higher energy. In contrast, the chain in Fig. 1(c) does 
not exhibit a clear peak at +Atip and also differs in the 
remaining subgap structure, with a resonance at around 
1.52 mV. Figure 1(e) and (f) show data for chains with¬ 
out a protrusion at their ends. These appear rarely 
and we can only provide data with reduced resolution 
330 /iV) due to inferior tip preparation. These chains 
also lack an unambiguous signature of a zero-energy res¬ 
onance, perhaps masked by the larger resonance which 
contributes considerable spectral intensity at Atip. 

Figure 1(a) shows that the peaks at opposite biases 
TAtip differ substantially in intensity. The same is ob¬ 
served in all other chains [40]. This is in contrast to 
expectations for Majorana peaks, which should be sym¬ 
metric, as a Majorana state has identical electron and 
hole wavefunctions [41]. This indicates that these peaks 
originate at least partially from trivial subgap states near 
zero energy. For a more detailed analysis, we focus on the 
chains shown in Fig. l(a,b). 

It is interesting to contrast the spectra for chains with 
those for individual Fe adatoms and dimers. Figure 2(a) 
shows the d//dV-signal of two species of single adatoms 
with different apparent heights. For both types, we ob¬ 
serve a single shallow Shiba state with energies 1.1 meV 
(type 1) and 1.2 meV (type 2), respectively. In con- 
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Figure 2. (a) d//dV-spectra of single Fe adatoms (two types: 

1 and 2). Both show a single, loosely bound Shiba-resonance 
close to the gap edge. They differ in apparent height by 
~ 20 pm at 50 mV, 50 pA. (b) Dimers show a variety of 
Shiba states with the lowest energy resonance at 150 /xeV. 
Scale bar corresponds to 1 nm. Atip = 1.39 mV. (c) d//dV- 
spectrum at the end of the chain (blue), compared to bare 
Pb(llO) (gray). The chain is the same as shown in Fig. 1(a). 
(d) Zoom on the spectrum in (c). The subgap spectrum shows 
a manifold of Shiba resonances. A peak at +Atip may be the 
fingerprint of a Majorana bound state. Low energy states lie 
at £o — 80 /iV (black solid arrow) and ei ~ 270 /xV (blue 
dashed arrow). The energies of the expected thermal reso¬ 
nances of low intensity are marked by their gray equivalents 
below Atip = 1.36 mV. 


trust, the dimer shown in Fig. 2(b) exhibits a richer 
subgap structure with a series of resonances reaching 
as low as 150 /xeV. We find that the subgap spec¬ 
trum varies in detail between different dimers, depending 
on interatomic distance, angle, and adsorption site [40]. 
This demonstrates strong coupling of Shiba states which 
can ultimately lead to the formation of Shiba bands in 
adatom chains. 

Figure 2(c) provides the data of Fig. 1(a) over a wider 
voltage range, with a zoom-in on the voltage range near 
Atip shown in Fig. 2(d). In addition to the peak at +Atip 
and a faint shoulder at —Atip, there is a nearby subgap 
peak at eo — 80 /iV (black solid arrows). These combine 
into a plateau-like structure near and just above TAtip. 
The data are also consistent with corresponding thermal 
resonances at ± (Atip — 5o) (gray arrows). The superpo¬ 
sition with a low-energy subgap resonance may explain 
the asymmetric peaks in Fig. 1(a). While a Majorana 
peak must be symmetric, conventional subgap resonances 
can be asymmetric, reflecting the asymmetry between its 
electron and hole wavefunctions. 
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Figure 3. Spatially resolved dJ/dlZ-spectra along chain (b) 
in Fig. 1, going from the onset of the Fe cluster (#1) to 
the bare surface (^35). A manifold of Shiba bands domi¬ 
nates the subgap structure. At the chain end (#29, same as 
Fig. 1(b)) and the Fe cluster (#1) a zero-energy resonance 
around ibAtip is visible. The lowest non-zero energy reso¬ 
nance has eo ~ 80 fieY in (#29), and is modulated along 
the chain. As a guide to the eye low-energy resonances be¬ 
low a bias of ± (Atip + 100 /iV) are marked by small ticks. 
The higher energy peaks at ~ 270 /xV shift in energy along 
the chain. Offset for clarity: —30 nS/spectrum. Distance 
between spectra: 0.33 nm. Atip = 1.42 mV. 


Further subgap peaks occur at higher energies, the 
next higher one at Si 270 /iV (blue dashed arrows 
in Fig. 2(d)). Similar peaks were identified in Ref. [18] 
as the coherence peaks of the induced topological gap 
(estimated at 200 — 300 /xeV). In view of the lower- 
energy peaks, this interpretation seems implausible for 
our chains. Instead, a Majorana-based interpretation 
would suggest that the lowest non-zero energy peak 
originates from the topological gap or is shifted above 
the topological gap by size quantization. This suggests 
that the topological gap is comparable or smaller than 

80 /ieV. 

While the Majorana states should be localized at the 
chain end, the topological gap is a property of the bulk 
spectrum and should be observable throughout the entire 
chain. Figure 3 shows spatially resolved d//dV-traces 
of chain (b) in Fig. 1. At the end of the chain (#29) 
[and the onset of the Fe cluster (#1)], there are peaks 
at TAtip as well as low-energy subgap resonances with 
So — 80 /ieV, similar to the chain shown in Fig. 2. The 
peaks at TAtip are again asymmetric. But while the 
zero-energy resonance is observed only at the end of the 
chain, the low-energy resonance at Sq is indeed observ¬ 


able elsewhere along the chain as indicated by ticks in 
Fig. 3. Notice, however, that the apparent energy of the 
resonance and its intensity varies widely along the chain. 

Many of the spectra are dominated by a strong sub¬ 
gap resonance at a bias of 2.3 mV. A false-color plot 
of the same d//dV-spectra reveals that its intensity os¬ 
cillates with a period of 2 nm and shifts slightly to 
lower energy in the center of the chain [see Fig. 4(a)]. 
The lower-energy resonances near Atip do not show such 
a clear periodicity. Interestingly, the variations appear 
correlated with the topography of the chain, as shown 
next to the false-color plot in Fig. 4(b). Apparent height 
and width of the chain show variations with a similar 
period of 2 nm, in agreement with Ref. [18]. 

Further insight can be gained from a false-color plot of 
d//dV-spectra along the chain in a larger bias range be¬ 
tween ±1.1 V [see Fig. 4(c)]. In the interior of the chain 
we observe two prominent features: a narrow resonance 
a around —380 mV and a broader resonance a' around 
450 mV. Similar resonances are present for all chains 
shown in Fig. 1, and are in agreement with Ref. [18]. 
These resonances decrease in intensity and finally disap¬ 
pear at the end of the chain, or close to the Fe cluster, 
respectively (see additional traces in [40]). We interpret 
these resonances as van-Hove singularities of the d-bands 
of the Fe chain. The simultaneous disappearance of a 
and a' suggests that these are the upper and lower edge 
of the same band crossing the Fermi level. Interestingly, 
resonance a shifts with the above observed periodicity of 
about 2 nm [40]. 

To understand the origin of the variations in the d- 
bands and the peaks in the subgap spectra, we have 
performed model calculations for a chain of Anderson 
impurities coupled to an 5-wave superconductor with 
spin-orbit coupling, based on Ref. [34]. We account for 
the modulations by including a potential which varies 
along the chain and reflects the local environment of 
the adatoms. We choose parameters such that one band 
crosses the Fermi level with band edges corresponding to 
a and a' and we assume strong adatom-substrate cou¬ 
pling. Following [34], we calculate the subgap local den¬ 
sity of states from a mean-field treatment of the impu¬ 
rity chain (see [40]). Figure 5(c) shows the subgap band 
structure of an infinite chain. It comprises two fiat parts, 
which can be viewed as Shiba states in the superconduct¬ 
ing substrate around ka 0 and ka tt, and a V-shaped 
dip that originates from the strongly dispersing impurity 
states. For direct comparison with the data, we com¬ 
pute the differential conductance. As temperature ex¬ 
ceeds the typical energy separation between subgap lev¬ 
els, we assume efficient quasiparticle relaxation which re¬ 
sults in dominant single-particle tunneling. In addition, 
we model the experimental resolution by introducing a 
broadening in the tip density of states. Note that these 
conditions preclude the observation of a quantized Majo¬ 
rana peak height [41]. Figure 5(a) shows the differential 
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Figure 4. False-color plot of the dJ/dV-spectra of (a) the subgap (same spectra as in Fig. 3) and (c) the d-band structure, 
aligned with the topography in (b). The subgap structure (a) exhibits variations of the Shiba peak energies and intensities, 
as well as modulations of the low-energy resonance at eo ~ 80 /xeV along the chain. A zero-energy resonance is found at both 
chain terminations. Resonances linked to the d-band structure (c) vary around —380 mV (a), and 450 mV (a') along the chain. 



Figure 5. Numerical results for a chain of 30 sites: (a) 
Color plot of the differential conductance along the chain at 
subgap energies for a superconducting tip. (b) Spatially vary¬ 
ing on-site energies of impurity levels, (c) Bulk subgap band 
structure of an impurity chain on the surface of an s-wave 
superconductor. We used Agampie/Atip = 0.958. See [40] for 
other parameters. 


conductance of a finite chain for subgap energies, includ¬ 
ing a spatially varying potential [shown in (b)] to model 
the variation of the d-level resonances in Fig. 4(c). The 
numerical results are consistent with key features of the 
experimental data: (i) The Majorana bound state has a 
short decay length of a few lattice sites as a consequence 
of the strong chain-substrate coupling [34]. (ii) Promi¬ 
nent peaks at eV = ±1.5Atip and ±1.9Atip signal the 
van-Hove singularities of the Shiba band. Their intensity 
modulations are correlated with the potential landscape 
of the impurity atoms. Here the effect of the corrugation 
is most visible, as the Shiba energy explicitly depends 
on the energy of the impurity level, (iii) The induced 
gap varies along the chain on atomic scales but is uncor¬ 
related with the potential landscape. Indeed, at strong 
coupling the induced gap only depends on the substrate 


gap and spin-orbit interaction and is insensitive to de¬ 
tails of the impurities. The fluctuations originate from 
finite-size quantization which is most visible at low ener¬ 
gies due to the low density of states in the V-shaped dip 
of the band structure [Fig. 5(c)]. 


Motivated by Ref. [18], we investigated the subgap 
spectra and possible Majorana signatures of Fe chains 
on a superconducting Pb(llO) substrate by scanning tun¬ 
neling spectroscopy. Using superconducting tips, a Ma¬ 
jorana state is expected to appear as a pair of resonances 
at TAtip with symmetric intensities. We associate the 
absence of this symmetry in the data with a nearby low- 
energy subgap resonance at 80 /ieV. Within a Majorana 
framework, it is natural to interpret this additional res¬ 
onance as the coherence peak of the induced topological 
gap which would then be significantly smaller than pre¬ 
viously believed. We show by model calculations that 
such an interpretation is in principle consistent with our 
observations. However, a conclusive confirmation of Ma¬ 
jorana end states in chains of magnetic adatoms would be 
greatly facilitated by experiments at considerably lower 
temperatures. Using superconducting tips at tempera¬ 
tures well below the induced gap might even provide ac¬ 
cess to the elusive conductance quantization of Majorana 
states [41]. 
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SUPPLEMENTARY MATERIAL 

I. d//dy-SPECTRA OF THE SUBGAP STRUCTURE OF SIX DIFFERENT CHAINS 

In Fig. l(a-f) of the main text we show d//dU-spectra at positive bias voltage of the subgap structure at the end of 
six different chains. For the sake of completeness, we show here the same spectra, but both positive and negative bias 
side [Fig. 6(a-f)]. All spectra exhibit an asymmetric intensity of the spectral weight at ±Atip, z.e., at zero-energy. 


T—'—I—•—I—•—I—•—n p—•—I—'—I—•—I—'—n p—'—I—•—I—'—I—•—n p—•—i—•—i—•—i—•—n p—•—i—•—i—•—i—•—n p—•—i—•—i—•—i—•—r 



- 2-10 1 2 - 2-10 1 2 - 2-10 1 2 - 2-10 1 2 - 2-10 1 2 - 2-10 1 2 


sample bias (mV) 


Figure 6. d//dF-spectra recorded at the end of six different chains. Same spectra as in Fig. l(a-f) of the main text, but here 
we present both positive and negative bias side. Chains in (a-d) are terminated by a small cluster, in (e,f) they have a sharp 
cut-off. The energy resolution in (e,f) is reduced (width of BCS resonance: ~ 330 /iV), due to non bulk-like superconductivity 
of the tip. Atip in mV: (a) 1.36, (b) 1.42, (c,d) 1.38, (e,f) 1.24. Chain lengths measured between chain end and cluster onset 
in nm: (a) 13.9, (b) 9.5, (c) 6.2, (d) 6.0, (e) 7.7, (f) 4.0. Setpoint: 5 mV, 250 pA. Lock-in: 15 912 Hz. 


H. Fe MONOMERS AND DIMERS 

In the main text we show d//dV-spectra of the two species of Fe adatoms, which are present on the (110) surface 
after the evaporation of Fe onto the cold sample {Tgampie < 10 A). Both species are adsorbed in between the [110] 
corrugation lines of the surface, yet they have different apparent heights {Az 20 pm at Vbias = 50 mV, / = 50 pA). 
In d//dV spectroscopy between ±3.5 V, we did not detect any d-state resonances. We occasionally found - besides 
single adatoms - also larger protrusions of different size and shape. As diffusion is hindered at the temperature of 
deposition, it is reasonable to assume that the vast majority of these are Fe dimers. We show the d//dV-spectra of 
three distinctly different dimers in Fig. 7. All of them show a rich subgap structure, indicating a strong interaction of 
the Shiba states. However, the subgap structure varies strongly, and depends on the inter-atomic distance, the angle, 
and the adsorption site of the Fe atoms. 
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Figure 7. (a-c) d//dV-spectra of various Fe dimers [reference spectrum on pristine Pb(llO) in grey as guide to the eye]. The 

spectra show diverse subgap structures, probably due to different adsorption configurations, (a) is the dimer shown in Fig. 2(b) 
of the main text. Atip in mV: (a) 1.39, (b) 1.38, (c) 1.37. Setpoint: 5 mV, 250 pA. Lock-in: 15 ^12 Hz. 
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III. d//dy-SPECTRA OF THE IRON CHAIN D-BANDS 

Figure 8 shows dJ/d^-spectra acquired on the chain presented in Fig. 4 of the main text. The spectra visualize 
the variations in intensity and energy of the resonances a and a' along the chain. At the protrusion at the chain 
end, these resonances have not yet developed. Instead, we observe a resonance around —930 mV (spectrum #40). It 
decays quickly along the chain. In spectrum (#38) it is hardly visibly anymore, while resonances a and a' start to 
gain intensity. Both increase in intensity when moving towards the center of the chain. In the center (^#20), a' is 
resolved as double-peak structure. This correlates with the tight-binding calculations presented in Ref. [18], and may 
be a hint that the broad resonance a' actually consist out of two resonances. The simultaneous appearance of a and 
a' suggests that they originate from the same band. Different positions within the chain show a shift of the resonance 
a by up to 150 meV, which is plotted in the bottom right of Fig. 8. At the Fe cluster both resonances decrease in 
intensity again (#1). 



sample bias (V) 


a (mV) 


Figure 8. d//dV-spectra (left) from top to bottom recorded along the chain (top right) starting from the iron cluster and 
going towards the chain end. The spectra are numbered from #1 to #50. The data is the same as plotted in Fig. 4(c) of the 
main text. Resonances a and a are marked by dashed lines. The energetic position of a is plotted versus the spectrum number 
and the distance (bottom right). The values are determined from a fit with a Gaussian peak and a linear background. Spectra 
are smoothed to the adjacent average of 25 points. Point distance 230 pm. Offset for clarity: —0.05 nS/spectrum. Setpoint: 
2 V, 850 pA. Lock-in: 2 mVrms, 912 Hz. 


A similar behavior of the state a is found in a chain, which is not terminated with an Fe cluster. As an example, 
we show a set of spectra in Fig. 9. The resonances a and a' again appear simultaneously close to the end of the chain. 
We do not observe a resonance at larger negative bias at the end of the chain, which resembles the one at —930 mV 
in Fig. 8. This supports the correlation of its appearance with the protrusion at the end of the chain. Though, 
approaching the large Fe cluster, a and a' vanish, and resonances at —800 mV and around 150 mV appear in the 
bottom of Fig. 9(e). Within the chain, resonance a exhibits a clear oscillation in its peak position with a periodicity 
2 nm, similar to the previously described chain. Resonance a' also shows a variation in its peak position, but not 
with a clearly identifiable periodicity. We also present the corresponding subgap structure in Fig. 9(a). The most 
prominent Shiba state at 2.17 mV shows an intensity oscillation, in accordance with the variation of a. Due to the 
limited energy resolution in these spectra, we cannot unambiguously conclude on an energy variation. 
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Figure 9. d//dy-intensity of the subgap structure (a), and the large bias range (c) as color plot with respect to the location 
along the chain and the sample bias. The topography in (b) is aligned to the position of the spectra in the color plot. The 
chain is the same as shown in Fig. 1(e) of the main text. The tip has a reduced energy resolution of 330 /iV. The position of 
a is plotted in (d). The values are determined from a fit with a Gaussian peak and a linear background, (e) shows the full 
dJ/dV-spectra of the data shown in (c). Data is smoothed to the adjacent average of 15 points. Point distance 3 A. Offset for 
clarity: —0.1 nS/spectrum. Setpoint: 2 V, 850 pA for (c) and (e), and 5 mV, 250 pA for (a). Lock-in: 2 mVrms, 912 Hz for 
(c) and (e), and 15 [N 912 Hz for (a). 


IV. DETERMINATION OF THE TIP GAP 

We use superconducting tips in order to improve the energy resolution beyond the Fermi-Dirac limit at 1.1 K 
and to detect asymmetries in the electron and hole components of the subgap states. d//dV-spectra show thus the 
measured spectral intensity of the sample convolved with the BCS-like density of states of the tip. A consequence of a 
superconducting tip with gap Atip is the shift of a sample resonance with an energy ±5 to a bias value of ± (Atip ± e) /e. 
The exact determination of “real” energies of the sample resonances thus relies on the correct determination of the 
superconducting gap of the tip Atip. 

Pb is a two-band superconductor with two gap parameters (Ai 1.42 meV and A 2 ^ 1.27 meV). They originate 
from two separated Fermi surfaces, and give rise to the double-peak structure in the d//dV-spectra [38]. The tip is 
prepared by controlled indentation into the clean Pb surface with a high voltage applied to the tip. This creates an 
amorphous superconducting Pb layer on the tip and yields a single gap parameter Atip, which is averaged over all 
directions. Depending on the layer thickness and quality, Atip can be similar or smaller than the bulk gap values. 

We cannot determine the parameters Ai, A 2 , and Atip independently from the BCS resonances in the spectra 
of pristine Pb(llO) alone. An independent determination of the full set of parameters (Ai, A 2 , and Atip) is only 
possible using spectra with a pronounced low-energy Shiba state, which gives rise to well-resolved thermal resonances. 
Fe dimers show such low-energy Shiba resonances (Fig. 7). The Shiba resonance and its thermal counterpart occur 
symmetric to Atip at ± (Atip + e) and ± (Atip — 5 ), respectively. This allows us to determine Atip unambiguously. 
Spectra of the pristine surface acquired with the same tip show clear BCS resonances at Atip + Ai^ 2 - Because Ai and 
A 2 are bulk properties of the substrate, their energy can then serve to determine Atip for every tip. This procedure 
enables a reliable determination of the energies of subgap resonances in each Fe chain. 
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(a) (b) 




Figure 10. Numerical results for the differential conductance measured with a superconducting tip at subgap energies along a 
chain of 30 sites, (a) Sample with a larger p-wave induced gap with kpa = 4.37r, same as the color scale plot in Fig. 5(a) of the 
main text, (b) Sample with a smaller p-wave induced gap with kpci = 8.37r, only the last 5 traces around Atip are shown. The 
dashed lines indicate eV = ±(Atip + Ap), with Ap the induced p-wave gap for an infinite chain. In (a) Ap = 0.096Atip and in 
(b) Ap = 0.045Atip. The remaining model parameters and the process of evaluating the conductance are the same. 


V. THEORETICAL MODEL 

We use the model introduced in Ref. [34], where the Fe adatoms are modeled as a chain of Anderson impurities 
hybridizing with the substrate BCS superconductor. The only difference is that we account for the corrugation of the 
chain by assuming a spatially varying on-site energy of the Anderson levels. We choose parameters such that only 
one spin band, say spin down, crosses the Fermi level, whereas the spin-up band is far below the Fermi level. As 
described in Ref. [34] we can evaluate the Green function within mean-field theory, from which we numerically obtain 
the local density of states (LDOS) at subgap energies. In our numerical calculations, the average on-site energy of 
the spin-down levels is = 200Asampie5 with roughly 10% spatial modulation (see Fig. 5(b) of the main text). The 
onsite energy of the spin-up states is = E^^i — 40000Agampie- We choose the bandwidth of the spin-down band 
as lOOOAsampie: the hybridization strength between the adatoms and the substrate superconductor F = 256Agampie, 
and the ratio Agampie/Atip = 0.958. Finally, we set the Fermi wavevector kpct = 4.37r and the Rashba spin-orbit 
wavevector kha = 0.267r with a the lattice constant. 

We start from the expression for the tunneling current (see supplement of Ref. [41]), 

^ ^ y* Tr (cj_) - G^’®®(a;)5f^(cj_)-h ^^(cj+)G^’^^(cj) - (cj+)G^’^^(cj)} , (SI) 

where R and L label the sample and the superconducting tip and all Green’s functions are located at the tunneling 
position. The electron and hole blocks are denoted by ee and hh. In the experiment, the tunneling current is 
dominated by single-particle tunneling events (cf. Ref. [38]) and we can approximate 

-iGU^) = A{u;)f{u;), (S2) 

iG^{uj) = A{uj){l - /(cj)), (S3) 

where f{uj) is the quasi-equilibrium distribution of the steady state and A{uj) = —2ImG^(cj) the LDOS of the sample. 
The distribution f{oj) satisfies the condition = 1 — /(cj), as in the case of the Fermi distribution function. For 

our numerical calculations, we assume that the system remains close to thermal equilibrium, f{oo) = npioj). Using 
the relations for the tip Green function g^{oo) = 27rip{uj)nF{oj), — —‘27rip{uj)nF{oo), we can write the current 

as 

I =—^ I dujp{u}-)Tr Ae{u})[nF{uj-) - f{u})] - p{uj+)T:r Ah{uj)[nF{oj+) - fF)] (S4) 

C 

=—fF J duj p{uj-)Tr Ae{co) [nF(w_) -/(w)], 


(S5) 
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Figure 11. Numerical results for the differential conductance measured with a superconducting tip at subgap energies along a 
longer chain of 200 sites with the same parameters as in Fig. 10(a), in particular the same modulation in the on-site energy of 
impurities, (a) Color scale plot of the conductance, (b) Waterfall plot for 30 sites in the middle of the chain. The dashed lines 
indicate the bias ±(Atip + Ap) corresponding to the p-wave gap of an inhnite chain. 


where denotes the 2x2 electron (hole) block of A. In the last line, we have used Tt Ae{uj) = Tt Ah{—uj) due 
to particle-hole symmetry, and the property = 1 — The conductance is then obtained by taking the 

derivative with respect to the voltage, 

27re^f^ C 

G{eV) = - - — J dco Tr{Ae(w)} {/o'(w-) [npiuj-) - /(w)] + p{uj-)np(L 0 -)} , (S6) 

where p'{ui) = dp{u)/du, and n'p{u) = dnF{ijj)/du. Inverting the sign of the bias voltage, we find 

G{-eV) = - — j dco Tr{A,{-ao)} {p'{ao.) [n^(a;_) - f{co)] + p(a;_)n^(a;_)} . (S7) 

Since in general Ti Ae{oj) ^ Ti Ae{—uj) = TTAh{uj), the conductance is not necessarily symmetric under reversal of 
bias voltage. However, an isolated zero-energy resonance yields symmetric peaks at eV = TAtip [41]. In Fig. 10(a) 
we show a waterfall plot of the differential conductance along the chain for the same parameters as the color plot 
in Fig. 5(a) of the main text. The dashed line indicates the p-wave gap for an infinite chain. The lowest-energy 
resonance in the finite chain is somewhat higher and varies along the chain due to finite size effects. In Fig. 10(b) 
we show that a smaller p-wave gap leads to an asymmetry at eV = TAtip. In Fig. 11 we show the conductance plot 
for a longer chain, where the lowest-energy resonance is now more homogenous along the chains and closer to the 
induced gap. The small offset from eV = ±(Atip + Ap) originates from the energy-dependent density of states in the 
superconducting tip and is not due to finite-size effects. 




















































